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Summary. Precise radial velocity measurements (Sv/c ~ 10~ 7 ) of Fell lines in 
damped Lya systems from very high quality VLT/UVES spectra of quasars HE 
0515-4414 and Q 1101-264 are used to probe cosmological time dependence of the 
fine structure constant, a. It is found that between two redshifts z\ — 1.15 and 
22 = 1.84 the value of Aa/ot changes at the level of a few ppm: (a Z2 — a Zl )/ceo = 
5.43 ± 2.52 ppm. Variations of a can be considered as one of the most reliable 
method to constrain the dark energy equation of state and improvements on the 
accuracy of the wavelength calibration of QSO spectra are of great importance. 



1 Introduction 

The late-time acceleration of the universe discovered from the luminosity 
distance measurements of high-redshift SNe la is now regarded as evidence 
for the existence of dark energy. In many models, the cosmological evolution 
of dark energy is accompanied by variations in coupling constants as, e.g., the 
fine-structure constant a [1]. Since theories predict different behavior of dark 
energy, from slow-rolling to oscillating, to study its dynamics the coupling 
constants must be measured with highest accuracy at each redshift. 

In order to fulfill this requirement we developed a method called 1 Single 
Ion Differential a Measurement, SID AM' [2, 3, 4, 5]. Based on the measure- 
ments of the relative line position of only one element Fell, it is free from 
many systcmatics which affect other methods, e.g., [6, 7]. It has been shown 
that SIDAM can provide a sub-ppm precision at a single redshift and that 
this level of accuracy is mainly caused by systematic errors inherited from 
the uncertainties of the wavelength scale calibration [5] (ppm stands for parts 
per million, 1CT 6 ). 
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Fig. 1. Histograms are the combined Fell profiles from the z = f.84 damped Lya 
system towards Q 110f-264 [8]. The zero radial velocity is fixed at z a b s = f.8389ff. 
The synthetic profiles are over-plotted by the smooth curves. The normalized resid- 
uals, (Tt al — J r ° bs )/cn, are shown by dots (horizontal dotted lines restrict the la 
errors). The dotted vertical lines mark positions of the sub-components. Bold hor- 
izontal lines mark pixels used to minimize \ 2 - The ranges at v < —50 km s _1 and 
at v ~ —30 km s _1 in the Fell A2600 profile are blended with weak telluric lines. 



In this presentation we consider our recent results [8] with the SIDAM 
procedure obtained from the analysis of very high resolution (FWHM ~ 3.8 
km s _1 , slit width 0.5 arcsec) and high signal-to-noise (S/N = 100-120 per 
pixel) spectra of Q 1101-264 (z cm = 2.15, V = 16.0). The observations were 
performed at the VLT Kueyen telescope on 21-23 February, 2006 under the 
ESO programme No. 076.A-0463. The total exposure time was 15.4 hours. 



2 Results 



The spectroscopic measurability of Aa/a = (a z — ao)/a at redshift z (ao 
refers to the z = value) is based on the fact that the energy of each line 
transition depends individually on a change in a [9] . It means that the relative 
change of the frequency u>o due to varying a is proportional to the so-called 
sensitivity coefficient Q = q/uio [4]. The q- values for the resonance UV tran- 
sitions in Fen are taken from [10], and their rest frame wavelengths — from 
[8,11]- 

The value of Aa/a can be measured from the relative radial velocity shifts 
between lines with different sensitivity coefficients. In linear approximation 
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Fig. 2. x 2 as a function of the velocity difference Av between the Fell A1608 and 
A2382/2600 lines [8]. The number of data points M = 305, the number of degrees 
of freedom v — 257. The minimum of the curve at Av = —0.18 km s _1 gives the 
most probable value of Aa/a = 5.36 ppm. The la confidence level is determined 



by A\ = 1 (dashed line) which g 1VGS (J Av = 0.08 km s or a Ac 



2.38 ppm. 



(\Aa/a\ < 1) we have [5]: 

Aa _ (f2 - vi) Av 

a ~ 2c(fii - Q 2 ) ~ 2cAQ ' [ ' 

where index T' is assigned to the line A 1608, and index '2' marks one of the 
other Fen lines (A2382 or A2600). 

The resulting normalized, vacuum-bary centric, and co-added spectra are 
shown in Fig. 1. Two independent data reduction procedures (ID and 2D) re- 
sulted in almost identical co-added spectra. We found a perfect fit of the Fe II 
profiles to a 16-componcnt model (shown by the smooth curves in Fig. 1): 
the normalized % 2 per degree of freedom equals \t = 0.901 {y = 257). The 
computational procedure was the same as in [5]. Since the Q values for the 
A2382 and A2600 lines are equal, their relative velocity shift, Ziw2600-2382, 
characterizes the goodness of wavelength calibrations, cr sca ic, of the corre- 
sponding echelle orders. We found that the value of Ziw2600-2382 = 20 m s _1 
is comparable with er sca i < 20 m s _1 estimated from the ThAr lines. So, in 
what follows we consider the A2382 and A2600 lines as having the same radial 
velocity, and calculate Av between this velocity and that of the line A1608. 

The radial velocity shift between the A1608 and A2382/2600 lines was 
derived by comparing synthetic profiles with their observed profiles and min- 
imizing x 2 ■ To find the most probable value of Av, we fit the absorption 
lines with a fixed Av, changing Av in the interval from —270 m s _1 to —90 
m s _1 in steps of 10 m s _1 (see Fig. 2). For each Av, the strengths of the sub- 
components, their broadening parameters and relative velocity positions were 
allowed to vary in order to optimize the fit and thus minimize x 2 (Av). The 
most probable value of Av corresponds to the minimum of the curve x 2 (Av). 
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In our case it is —180 m s _1 . The lcr confidence interval to this value is given 
by the condition Ax 2 = X 2 — Xmin = 1 (^ nc horizontal dashed line in Fig. 2. 
It is seen from the figure that a a v — 80 m s -1 , or Aa/a — 5.36 ± 2.38 ppm. 

We have specifically investigated those systematic effects which could in- 
troduce a non-zero difference between the blue and the red lines (A1608 and 
A2382/2600, respectively) and thus simulate a variation of Aa/a at the ppm 
level: (1) possible isotopic shifts caused by unknown isotope abundances, (2) 
effects of the unresolved components, and (3) possible blends. We also com- 
pared our previous results [8], where Fen lines falling in both UVES arms 
were used, with the present measurements to check up possible shifts caused 
by different changes of isophote onto the slit during integration and did not 
reveal any of them. We found that the largest systematic error does not 
exceed 30 m s _1 , i.e. (JAa/a.sys < 0.89 ppm [8]. 

We note that the accuracy of (JAa/a — 2.38 ppm represents a factor of 

I. 5 improvement with respect to our previous result (JAa/a = 3.8 ppm [4] 
obtained from the archive data of lower resolution (FWHM ~ 6 km s _1 ) 
but the same S/N ratio. Thus the higher spectral resolution is shown to 
significantly contribute to higher accuracy in the Aa/a measurements. 

The comparison of Aa/a = —0.07 ± 0.84 ppm at z\ = 1.15 towards 
HE 0515-4414 [5] with the measured quantity Aa/a = 5.36 ± 2.38 ppm at 
z 2 = 1.84 shows a tentative change of the value of Aa/a between these two 
rcdshifts: (a Z2 — a Zl )/a a = 5.43 ± 2.52 ppm. 
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